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1. INTRODUCTION 

There has recently been a great deal of interest in QCD at high temperature, in part 
due to the hope that a new state of matter, the quark gluon plasma, might soon be observed 
in relativistic heavy ion collisions [1]. There are, however, very few observable properties of 
this plasma that can be predicted directly from the fundamental theory. Most experimen- 
tal properties are derived from phenomenological models that are often ad hoc and contain 
many parameters. The few physical quantities that can be calculated using perturbative 
QCD have recently been the subject of much controversy [2], in part due to the apparent 
gauge fixing dependence of the gluon plasma damping constant at one loop[3]. This gauge 
fixing dependence has been attributed by some to the supposed breakdown of perturbative 
QCD [4], and by others to the need for a manifestly gauge covariant or invariant linear 
response theory [5]. The correct resolution to this puzzle has recently been proven to lie in 
the innacuracy of the naive loop expansion as a self-consistent approximation scheme: an 
accurate, and gauge independent damping constant can only be obtained by resumming an 
infinite number of loop diagrams. Although the need to re-sum was realized fairly early [6], 
an explicit prescription for performing the necessary resummation has only recently been 
provided by Pisarski[7]. The proof that the resummation does indeed yield a gauge inde- 
pendent damping rate was provided using two distinct methods. The first, due to Braaten 
and Pisarski (BP)[8], involved explicitly performing the resummation required and using 
Ward identitities satisfied by the improved perturbative n-point functions to show that the 
same resummed damping constant is obtained to leading order in strict Coulomb gauge 
and in all covariant gauges. A similar result was also obtained in covariant gauges using 
a partial resummation technique[9]. The second method, due to Kobes, Kunstatter and 
Rebhan (KKR)[10] used properties of the generating functional for QCD to provide a non- 
perturbative proof of the gauge fixing independence of the locations of the physical poles 
of the gluon propagator. 

In effect, both the BP and KKR proofs were based on the fact that the gauge depen- 
dence of the physical structure functions vanishes "on- shell" . That is, given the structure 
function f{k) whose zeroes determine the location of the physical poles in the gluon prop- 
agator, one can show that 

Afik) = fik)AXik) , (1.1) 

where AX{k) can be calculated perturbatively. Eq.(l.l) has the immediate consequence 
that the physical pole position (and hence also the mass and damping rates, which are 
essentially the real and imaginary parts of the pole) are gauge fixing independent, providing 
that AX{k) does not have poles when f{k) = 0. KKR provided a non-perturbative proof 



of (1.1) for gauges theories in general, and then discussed how to apply the result to a 
self-consistent approximation for the gluon plasma damping rate. BP on the other hand, 
provided a perturbative proof of (1.1) using the resummed Feynman rules appropriate for 
finite temperature QCD. Clearly, the former proof was very general and necessarily quite 
formal, while the latter proof considered only a restricted family of gauges. Moreover the 
proof of Braaten and Pisarski was also formal, to the extent that it did not explicitly 
calculate the particular /S.X{k) that arose, but instead assumed that it would be well 
behaved. KKR argued that the pole structure of /S.X{k) would be essentially determined 
by the poles of the full ghost propagator, and hence distinct in general from the solutions 
to the physical dispersion relations, but as will be shown below, their analysis must be 
re-evaluated in the presence of infrared divergences. 

The purpose of the present paper is to extend the Ward identity analysis of Braaten 
and Pisarski to display explicitly the gauge fixing dependence of all the components of the 
resummed gluon self energy and to include all possible linear gauges in which rotational 
symmetry (in the rest frame of the heat bath) is unbroken. Thus we consider all covariant 
and non- covariant gauges in which the four-velocity of the heat bath is the only additional 
vector that appears. 

The paper is organized as follows: In Section 2, we present the Feynman rules for 
one loop QCD calculations in the class of gauges considered. This serves to establish con- 
ventions and notation and also to remind the reader of the source of the tree level Ward 
identities obeyed by the bare n-point functions. Section 3 explains briefiy the need for con- 
sidering higher loop diagrams, and describes how a corrected damping rate can be obtained 
using the resummation technique proposed by Pisarski [7] and implemented by Braaten and 
Pisarski [8]. Section 4 uses the Ward identities to analyze in detail the gauge dependence 
of the gluon self energy, and its associated damping rate, in this improved approximation 
scheme. In Section 5 we examine the integrals that must be evaluated in Covariant and 
Coulomb gauges in order to prove gauge fixing independence of the damping constant 
and find the surprising result that in covariant gauges, the imaginary part of AX(/c) does 
develop poles on the physical mass shell and leading to gauge dependent damping rates 
unless an infrared regulator is maintained throughout the calculation. This extends to glu- 
ons a similar result concerning the gauge dependence of the resummed fermion damping 
rate[lll. Section 6 closes with conclusions. 



2. PRELIMINARIES 

For simplicity we consider only gluons. The discussion can be generalized to include 
quarks without much difficulty. The Feynman rules for the theory are derivable from a 
generating functional of the form: 

^ = / n dA';^{x)dC''{x)dC''{x) eMlci + Ig.f. + hh) ■ (2.1) 



a,iJ.,x 



In the above, Id is the classical SU{N) Yang-Mills action: 



Id = -7 / d'xF'^^-'F^, , (2.2) 



where the colour indices a, 6, c. run from 1 to A^ — 1 and the Yang-Mills field strength is 

F;, = d,A: - d^Al - g^^AlAl . (2.3) 

The classical action is invariant under the following infinitesmal gauge transformations of 
the vector potential: 

M«(x)=I)^^(x)A^(x), (2.4) 

where we have defined the covariant derivative operator: 

D';!'{x)^d-%-gr'^A';,{x) (2.5) 

and A"(x) is an arbitrary, infinitesmal Lie algebra valued field. The fields C and C" are 
the Grassman-valued Faddeev-Popov ghosts necessary to make the source-free generating 
functional independent of the gauge fixing term /g./., which is of the form: 



Ig.f. = --. I dM^F'^f , (2.6) 



1 
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where JF" is, in principle, an arbitrary function of the vector potential, restricted only by 
the condition that its associated Faddeev-Popov operator be invertible. We will restrict 
the consideration to linear gauges 

jra=jrM(^)^a(^) ^ (2.7) 

where JF^ is an arbitrary (possibly non-local) differential operator, in which case the 
Faddeev-Popov operator takes the form: 

Q"' = jA^fT^^^ = •^^(^X(^) ■ (2-8) 



The ghost action is then : 

Igh = f d'^xC"j'^D''J'C^ . (2.9) 

Since we are interested in finite temperature Green's functions [12] using the imaginary 
time formahsm, we wiU assume that we have aheady Wick rotated to Euchdean space-time, 
with metric signature (H — 1-++). AH functions are therefore periodic in imaginary time with 
period /3 = 1/T, as required by the Matsubara formahsm. This in turn requires the energy, 
/cq, to be an integer multiple of 27rT. The retarded real time thermal Green's functions 
from which the plasma parameters can be extracted are obtained from the imaginary time 
Green's functions by analytically continuing the external momenta back to real time, after 
all Euclidean sums and integrals over internal momenta have been performed. 

At finite temperature, the velocity of the heat bath introduces a preferred direction in 
space-time which breaks manifest Lorentz invariance. We will denote this velocity by the 
four vector n^, and assume that we are in the rest frame of the heat bath, so that n^ = ^o- 
If the gauge fixing condition does not break rotational invariance, then the momentum 
space representation for the most general JF^ is: 

J^^ = d{k)k^ + b{k)n^ (2.10) 

for arbitrary functions d{k) and b{k) of the four momentum k^. 

In order to simplify the following discussion it is useful to introduce some relevant 
tensors. First of all, we will need the projection operator: 

^; = ^; - ^ ^ (2-11) 

where /c^ = (k^)'^ + k'^. This can be used to isolate the component of n^ that is orthogonal 
to the four momentum, /c^, namely: 

h^{k) = Pl^rf = 5^- k^"^ . (2.12) 

Note that n^ = /c^//c^. It will be useful to re-express the momentum space gauge fixing 
condition in terms of the resulting orthogonal basis: 

T^ = a{k)k^' + h{k)fi^, (2.13) 

where a{k) = d{k) + h{k)n • kjk^ . 



Finally, we define: 

J^flU — p^.u _ "- "- 



•n? 



n A _ Mi ) ' (2-14) 

which projects out spatially transverse modes, and is orthogonal to both k^ and h^. 

The corresponding Feynman rules are as follows: 
The bare 2-point function is 

j^^(o)-i j a6 _ ^ab A 2p^^ ^ 1 (^-^ ^ afc^)(6n, + aK)\ , (2.15) 

with associated propagator: 

6'^Dfiik) = r' (Ip,. - -|j(s,fc. + s„*,.) + {^^^^) ^) . (2,16) 

The ghost propagator and ghost-ghost-gluon vertex in this family of gauges are, respec- 
tively: 

and 

gP^'^T^'ik) = gr^^iak^ + hh^") . (2.18) 

Due to the linearity of the gauge fixing condition, the three and four point functions depend 
only on the classical action, and are, respectively : 

-igr^^r^^f'ip, q, k) = -igr^^{{p - qY6^'' + perm) , (2.19) 

and 

-g^NT^^'P^ip, q, r, s) = -g^N{2d^''dP^ - d^Pd''^ - 5^^5^^) , (2.20) 

where p + q + k = and we have traced over the colour indices in the four point function. 
An important aspect of the above bare n-point functions is that they obey the following 
Ward identities (suppressing colour indices): 

k^r^'^ik) = , (2.21a) 

kpF^'-^ip, q, k) = -V^^{p) + V^^{q) , (2.216) 

s^T^'^^'^ip, q, r, s) = T^^^ip + s, q, r) - T^'^^ip, q + s,r) , (2.21c) 
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where T^'^ = k'^P'^'^ refers to the gauge independent part of the tree level 2-point function 
in (2.15) above. These identities follow directly from the gauge invariance of the classical 
action from which the n-point functions were derived. In particular, one can express the 
gauge invariance of the classical action in terms of the following functional differential 
equation: 

which expresses the fact that the directional derivatives of the action along the gauge 
orbits are zero. The Ward identities then can be derived by functionally differentiating this 
equation with respect to the vector potential, and then evaluating the resulting expression 
at A" = 0. Differentiating once yields (2.21a), twice yields (2.21b) and three times yields 
(2.21c). 



3. RESUMMED PLASMA PARAMETERS 

The loop expansion does not provide a self-consistent approximation for the finite 
temperature QCD damping rate. In general, an infinite number of loop orders contribute 
to the leading order result[6]. This can be understood by noting that at finite temperature, 
QCD contains two independent mass scales: T and gT «T. Momenta that are of order 
T or greater are called "hard" , while those that are of order gT or smaller are called "soft" . 
For n-point processes in which all legs are of order gT, one loop corrections in which all 
internal momenta are hard contribute to the process to the same order as tree level[7]. 
Moreover, these are the only such large "corrections" . This implies, for example, that the 
dispersion relations can be consistently calculated to order g'^T'^ by considering only one 
loop self energy diagrams with hard internal momenta. No higher loop corrections are 
required. Explicit calculations show that the hard loop self energy diagram is real. The 
damping constant therefore appears only at the next order in the resummed perturbative 
expansion, namely g^T. One loop calculations to this order require summing over diagrams 
with soft internal legs and as mentioned above, for soft processes (internal or external) one 
needs to encorporate the hard thermal loop corrections to the vertices and the propagator. 

It is worth stressing that this algorithm is only valid for the leading order (g'^T) 
damping constant, and must be modified in order to calculate corrections either to the 
dispersion relations, or the damping rate. For example, in order to calculate the real part 
of the dispersion relations to order g'^T, it is necessary to consider higher loop diagrams 
involving the modified n-point functions [7]. 

Remarkably, the hard thermal loop corrections to the n-point functions are gauge 
fixing independent and obey the same Ward identities as the tree-level n-point functions 
[7,8,13]. ^ Let dU^'^ , SF'^'^p, and ST'^'^p'^ denote the hard thermal loop corrections to the 
two, three and four-point functions, respectively. Then it can be shown that [8, 13]: 

k^*r^''{k) = , (3.1a) 

A;/r'^^^(p, q, k) = -*r^'^(p) + -^T^^iq) , (3.16) 

s,*r^'^^^(p, q, r, s) = *r^'^^(p + s, q, r) - *r^'^\p, q + s,r) , (3.1c) 

where we have defined the corrected n-point functions: 

*rM^ = k^P^'^ + OT'^^ , (3.2a) 



^ The gauge invariance of the hard thermal loop contributions to the quark self energy 
was first shown by Klimov and Weldon[14]. 



The tree level Ward identities obeyed by the *r imply that the resummed n-point 
functions can be derived from a gauge invariant (resummed) effective Lagrangian. The 
existence of this effective Lagrangian, which has explicitly been constructed in Refs. [13,15], 
plays an important role in understanding the connection between the BP and KKR proofs 
of gauge fixing independence. 

In terms of the *r n-point functions, the resummed expression for the leading order 
imaginary part of the real time gluon self energy is: 



Im *n'^"(Ko, 1^1) = ^Im [ dp *D^p *T^''''f'{k,-k,p,-p) 

^ Jsoft 

^ ^-Im [ dp *r«^^(p,Q,A;) *D^^,{p) "D^p^iq) *T"'f'''{p,q,k) 

Jsoft 



2 .'soft 



2 f {a{p)p^ + b{p)n^'){a{q)q'' + b{q)n''{q)) 

Jsoft a{p)p^a{q)q^ 

where *D^^ is the effective propagator defined as: 

^ {*V^i, + -{hh^ + ak^){hh^ + ak„))~ . (3.4) 

In the above and what follows, Im refers to the imaginary part of the integrals after 
analytic continuation of the external momenta: iko ^ Kq + ie, and 

dP^T E / 703' (3-5) 

-oft po=2-nT-^Wt i^TT) 

where the spatial integration is performed over soft momenta only, and therefore must be 
cut-off at some value A such that gT « A « T. In the calculation of the imaginary 
part of the self energy, which is all that we will be concerned with here, the precise value 
of the cut-off turns out to be irrelevant. 

Standard linear response theory [16] implies that the response of the plasma to small 
external perturbations is determined by the physical poles in the retarded thermal gluon 
propagator. It is convenient to write the full propagator in terms of the hard thermal loop 
propagator, *-D^;y, (suppressing colour indices) and the resummed self energy *Il^i^ (with 
the imaginary part at order g^T given in (3.3)): 
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i)^- = (*d;J+ *n^.)-' 



1 1 



T{k) 



C{k) 



1 1 k^k'^ 

n^'fi" + TTTVTin^'k'' + h^k^") + 



C{k) 



D{k) fc2 _ • 



(3.6) 



For the class of gauges given in Section 2, the most general form of the gluon self energy 
in momentum space is : 



^J.~u 



n^n 



U^" = *Ut{k)A^'' + *Ui{k)^^ + *ne(A;)(fi'^fc^ + n^k^) + *Ud{k) 



TV 



k^'k'' 



(3.7) 



Note that it is not transverse in general. The components of the propagator in (3.6) can 
be related to those of the self energy as follows[17]: 



T{k) = k'^ + *U 



t , 



c{k) = kl + *n, + 



2^4 






(l+^*nd/a2p)_ 



C{k) 
D{k) 
where we have defined 

and 



g jk^ + r'^d/a^) 
bn^ (1 + ^ *Ujab) ' 

a^k^ ( l + i *nd/A;2a2 



i \kl+ *Ui+fi%^/C 



kl = k'^ + SUq 



ki =k^ + 6Ui 



In terms of the above decomposition: 



*r^'^ = klA^" + ki 



fi-iy 



n'^n 



L ~ A 1 



(3.8a) 
(3.86) 

(3.8c) 
(3.8rf) 

(3.9) 
(3.10) 

(3.11) 



consistent with the Ward identity (3.1a). 

The explicit expressions for the hard thermal loop corrections to the two-point function 
are, after analytic continutation back to real time z/cq — ^ Kq + it [8,14]: 



OTT(i^o,|A;|)~4 



2P 



1 



Kl 2\k\ \Ko-\k\ 



(3.12) 
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and 



dUL{Ko,k^)r^+3m 



^ Kn Ko + \k\ 

1 ^In 



(3.13) 



m \Ko-\ 

where m'^ = g'^NT'^ /9 in the absence of quark fields. The hard thermal loop corrections to 
the three and four-point functions are not required in the following discussion. In fact, as 
shown by Braaten and Pisarski[8], and elaborated below, the specific form of the corrected 
n-point functions is irrelevant to the formal proof of gauge independence of the resummed 
plasma damping rates. It is necessary and sufficient that the n-point functions satisfy the 
Ward identities (3. la)- (3.1c). 

The physical modes of the gluon in the plasma are transverse with respect to the 
momentum four-vector. As can be seen above, the presence of the preferred four-vector 
n^ leads to the existence of two independent physical modes at finite temperature . The 
first, determined by the coefficient of A'^'^ is the usual spatially transverse gluon mode. 
The second is a spatially longitudinal collective mode which is absent at zero temperature 
and is usually called the"plasmon mode". The dispersion relations and damping rates for 
the gluon and plasmon modes are therefore obtained from the real and imaginary parts, 
respectively of the solutions to: 

T{k) = , (3.14a) 

C{k) = . (3.146) 

Note that in principle, the non-transverse parts of the self energy can contribute to 
the longitudinal dispersion relation and damping rate. However, in linear gauges, *Ilc and 
*Ild are not independent, but are related by the following (exact) Ward identity [18]: 

T^J^.D^" = C , (3.15) 

which, in terms of the parametrizations given above reduces to: 

*Ud{kl + *Ui) = n^e *n^ . (3.16) 

Eq.(3.16) can be used to express the longitudinal structure function in the following 
form [17]: 

C(k) = ^ f^, ^ . 3.17 

(kl + *n. + 5^ *n^) 

Thus, in general, it appears that *Ilc can not only shift the location of the longitudinal 
pole, but also change the order of the pole from first order to second order. One of the 
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results of the present paper, however, is to show that in hnear gauges, *nc does not in 
fact contribute to the leading order dispersion relations. This has previously been only 
verified in strict Coulomb and covariant gauges [8]. 

Given the two structure functions T and £, it is straightforward to derive the disper- 
sion relations K^ = uj^'\k), which are determined by the zeros of the structure functions. 
Normally one assumes that the system is underdamped, so that the imaginary part of the 
pole is small compared to the real part. That is: 

where 7 << ujp. In this case, given a general structure function of the form: Kq— f(Ko, /c), 
the dispersion relation becomes: 



ul = Ref{up,k) , (3.18) 



while the damping rate is: 



7 



Imf{uj, k) 



, ,2 9 



Ref{oJ,k) 



(3.19) 



J LO^=LOr. 



For the gluon excitation at rest {\k\ = 0), Eqs.(3.12) and (3.13) yield the same mass 
for the longitudinal and transverse modes, namely 

ml,L = ml . (3.20) 

Application of the above techniques in Coulomb gauge yield a damping rate of [8] : 

7T,L ~ 6.63 ^-— . (3.21) 

z47r 

In reference [8] the Ward identities were used to prove that the same damping rate would 
result in covariant gauges, while Ref.[10] provided general arguments concerning the gauge 
independence of 7. In the following section we examine explicitly the gauge dependence 
of the imaginary part of the leading order resummed self energy (3.3) in a general class of 
linear gauges. 
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4. GAUGE DEPENDENCE OF THE RESUMMED SELF ENERGY 

We will now calculate the gauge dependent part of the imaginary part of the resummed 
self energy as given by Eq.(3.3) above. For convenience we write the propagator (3.4) in 
the following form: 

where we have defined 



a{k) = \^ + k'/3' . (4.3) 



Since the modified 3- and 4-point functions are gauge fixing independent, all the gauge 
dependence in the self energy will come from the gauge dependent terms in the propagator. 
This can be isolated by defining: 

*]jt^'^ ^ *]j>niy j^ /\*]jt^'^ ^ (4.4) 

where we consider the gauge variation about an arbitrary "fiducial" gauge with resummed 
propagator: 

k^T kl kl' ^ kl 

Note that the above form of the propagator assumes that ^o = in the fiducial gauge, so 
that the gauge fixing condition is enforced by a delta function in the path integral. For 
example, if the fiducial gauge is strict Coulomb gauge, then (3o = ko/k'^, while in strict 
Covariant gauge, (3q = 0. These examples will be treated in detail below. The gauge 
dependent part or the propagator can now be written: 



where A(3 = (3 — (3q and 



v 



kl kl k^ 



Aa^k'{(3'-(3',) + %. (4.7) 
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The expression for the resummed self energy *n^'^ in Eq.(3.3) can easily be separated into 
its gauge independent (for a specific fiducial gauge) and gauge dependent parts. We con- 
sider each term in expression (3.3) in turn. The first term, involving the 4-point function, 
has the following gauge dependent contribution: 



Im < A*D *r(4) > = ^Im f dp A*L>a/3 *r^" 

2 J soft 

g'N 



ap 



-Im 



dp 



'soft 

A« PcPp 

-»2 9 

PI p 



2—{na{p)Pf3 + np{p)pc,) + 



Pl 
T^^^^ik^-k^p^-p) 



g^NIm [ dp ^^fiMi *r^""(A;,g,p) + *r"^"(A;,Q,p)) 

J soft P L 



g'N 



Im 



dp 



-»2 
PI 

Aa 



soft ' PlP'^ 



(2 *r^'^(A;) - *r^''{k-p) - *r^''{q)) 



(U 



where the last line was obtained by applying the Ward identities (3.1b) and (3.1c) succes- 
sively. Note that A/3(-p) = -A/3(p) and Aa{-p) = Aa{p). 

The gauge dependent contribution from the 3-point functions (2nd term in (3.3)) 
is considerably more complicated. It can, however, be simplified by applying the Ward 
identities whenever an internal 4-momentum is contracted with one leg of either a three, 
or two point function. The result can be summarized as follows: 



Im < *T^'^^ A*D *D *r*^"^'' > + < *r*^'^'' *£> A*D *t''^'^ > + < *r*^"^'' A*D A*D *r*^"^'' > 

+ g'^NIm [ dp ^^nf3{p) *r''^^{q,p,k) + {p ^ u) 
J soft I Pl 

Aa{p) 



+ g^NIm f dp -^P^ *T'"'{q) 

J soft P lP 



+ g'^NIm / dp 

J soft 

Im, I dp 

J soft 



- A(3{p)h^{p) (q- (3o{q)fi-{q) \ , , . 

A^{p)Am r~ 



v?{q) 
h^{p)7V^{q) + {fx ^ u)] 



2 ' .L^f/'' 'n?{p)ii?{q) 

where we have used the identity: 

h^{k) *r^"(A;) = fc|n7n2 
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(4.9) 



(4.10) 



and defined: 



J soft q L \ IP "■ ipj J/ 



oft qL 

-g^NImf dp^ 

J soft qiq 



p,.«„),,„„,^ 



Loft qW ^'^(p) 

g^N^ f ^ AP{p)Al3{q) 



Im [ dp ^^^fJt^^^^ (nAp)np'{q) *T-'^'^{p,q,k)p^ 
J soft P r.q L ^ 



2 J soft PlIl 

q'^N f 
+ ^-^—Im / dp AP{p)AP{q) 

2 Jsoft 



n^{p)nf^{q) n'^{q)n"{p) 
n^{p)ql n'^{q)pl 



(4.11) 



and 



Y^^, = g'NIm / dp ^^ *I^L'(p) 

soft q Lq 



Jsoft Piqi 

+ g'^NIm / dp — ^^ ^2 2 — [^a{p)Pc^' + f^a'{p)Pc 

Jsoft P tM tM 



g^N f Aa{p)Aa{q)pcPc 



+ ^ I^ / ip — ' -— - . (4.12) 

^ Jsoft p lp q Lq 

Finally, the gauge dependent part of the ghost contribution is 

Im < Aghost > = 2g^NIm f dp ^§^fi''{q) (^ - /3o(p)?^ 

Jsoft ^ [q) \p ^ [p) 

-g'NImf dp^^n^ip)n^q). (4.13) 

Jsoft ^ (p) n iq) 

It should be noted that in order to derive the above expressions, extensive use was made 
of the identity (4.10) and: 

Im / dpf{p,q) = Im / dpf{q,p) , (4-14) 

J soft J soft 

which is valid providing the cut-off is much greater than the external momentum. 

After summing (4.8), (4.9) and (4.13), one finds a complete cancellation of all terms 
that do not contain at least one factor of the gauge invariant two point function *T^^{k). 
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The final expression for tlie gauge dependent part of the imaginary part of the self energy, 
therefore takes the simple form: 

im A*u'"'{k) = z{k) *r^''{k) - *r^^x^ - x^ *r"^ - *r^"yaa' *r"'^ , (4.15) 

where we have defined 

Z{k) = -g^NIm [ dp ^^ , (4.16) 

J soft PIP 

and X^" and Y'^'^ are given in (4.11) and (4.12) above. It should be stressed that Eq.(4.15) 
above is valid for finite and/or infinitesmal gauge changes. 

Expression (4.15) contains the desired information about the gauge dependence of 
the resummed self energy. First of all, by contracting with k^k^, using the decomposition 
(3.7) and the transversality of *r^'^, one can easily verify that the gauge dependent part 
of *Ild is identically zero to this order. This is consistent with the Ward identity (3.16) 
which implies that, if *nc ~ g^T, then *Ild must vanish to this order. 

Secondly, by contracting Im A*Il^'^ in (4.15) with k^fii, one obtains: 

Im A*U, = kl^h^X^^k, . (4.17) 

To derive (4.17) one needs to use the transversality of *r^'^ , and the identity (4.10). Since 
Braaten and Pisarski verified that /m *Ilc vanishes on the longitudinal mass shell in strict 
Coulomb gauge, we have proven that this will be true in any gauge, providing only that 
the coefficient of /c^ in (4.17) does not have a pole at kj^ = 0. On this basis one can 
conclude that the non-transverse part of the resummed self energy will not contribute to 
the damping constant of the longitudinal mode in any gauge. 

By contracting (4.15) with the projector A^^, one finds that: 

Im A*Ut{k) = klZ{k) - klA^^X^'^ - i^A^^Y'^'' , (4.18) 

while contracting with n^fii, yields 

Im A*Ui{k) = klZ{k) - 2fci^^i^^ - (^i)'V^'-^^'^ . (4.19) 

Thus, the gauge variation of the structure functions vanishes on shell, providing that the 
corresponding coefficients in the above equations do not develop poles on the mass shell. As 
discussed in the introduction and in Ref.[10], this implies the gauge fixing independence 
of the resummed damping constants, and generalizes the formal proof of Braaten and 
Pisarski [8], who considered only Coulomb and Covariant gauges. In the next section we 
examine in more detail the integrals that appear in (4.18). 
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5. DAMPING RATE GAUGE DEPENDENCE FOR COVARIANT 
AND COULOMB GAUGES 

Let us study the gauge variation of the transverse structure function as an example. 
We will evaluate explicitly the coefficients of k'^ in (4.18) in order to check their behaviour 
on mass-shell {k^ = 0). 

Covariant Gauges 

We start by investigating variations around a fiducial covariant gauge, taking 

A/3 = (/3 = /3o = 0) , Aa = e^ 

(with a trivial change of notation for the gauge parameter as compared with (4.7)). 

Using the expressions for X^'^ and Y^'^ as given in (4.11) and (4.12) and the following 
identities: 

A,.(fc)P-(p) = 2-^ + i|^ 

(k-p)^ 



k^p ^ 
A,Ak) n'{p)n^{p) = ^ f 1 - ^' 

p y k^p ^ 

A,Ak)p^p''=p'-^-^: 



(5.1) 



where p"^ = Po + p'^^ ^^ get 



Im A*Ut = - ig^N Im 



.2 



dp I p^ {k ■ p)"^ 

soft ^ 



soft "^ \ -f- 



4 ^y ?v^ I fc2-.2 

soft ^ 



2 
2 



It should be noted that, in the derivation of (5.2) from (4.18), the analytic continuation of 
the external momentum is implicitly assumed so that Im A*nt(/c) = Im A*nt(i^07 k) . 
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In the following, we evaluate the two integrals in (5.2) that are potentially dangerous 
on-shell: Ii{k) and /2(^) denoted respectively as "transverse" and "longitudinal" and given 
by: 

1 FT 

soft 

J Q P Pl 

soft 

The integrals in (5.2) which do not depend on the plasma excitation properties pt 
and PL, respectively, are treated in Ref.[19]; indeed they are well behaved on the physical 
mass shell, kr = 0. 

(i) Calculation of the transverse integral 

We restrict ourselves to the calculation for the gluon excitation at rest: k = 0. In 
order to evaluate the double pole we use the identity 

— - = lim I — — — - 1 -r . (5.4) 

It is easy to derive that on the mass shell, i. e. when Kq = nig, the integral /i has no 
imaginary part except for 171 = 0. Therefore we define li{mg,0) by taking the following 
limits: 

/i(m„,0)= lim lim < ( -- — ^ ) /m / —^ ^r^ / ^ (5-5) 

K soft ) 

with p = {w,p), q = (w', —p), w' = Kq — w. 

Following the method of Refs.[8, 20], which we already used in a related study of the 
fermion damping rate [11], we calculate the imaginary part above by using the spectral 
representation of the effective transverse gluon propagator At{p) = -\r [14, 21] and of 
the free one A(g) = i}^i , as 



(g^ + m^^pi^ 

TT /" dp 



soft 

+ 00 

dw dw' 5{Ko — w — w')nB{w)nB{w')pt{w,p)e{w')5{w''^ — E'^] 



nB{Ko)J (27r 



)3 



-oo 



(5.6) 
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with E"^ = p^ + m? and pt{w,p), the spectral density associated to At(p) given in Ref. 
[21]. The spectral density associated to A(g) is e{w')5{w''^ — i?^), and ub is the Bose 
distribution. 

Integrating with respect to w' and using the symmetry property: 



pt{-w,p) = -pt{w,p), 



we get 



1 



CXD „^„2 pOO 



pdp 



dwpt{w,p) 



SirnBiKo) Jq E Jq 

■ {nB{E)[nB{w)6{Ko - w - E) + {1 + nB{w))5{Ko + w - E)] 

+ (1 + nB{E))[nB{w)5{Ko - w + E) + {1 + nB{w))5{KQ + w + E)]}. 

We first concentrate on the time like pole contribution in pt{w,p): 

pr''' = priw.p)s{w-wtip)), 



(5.7) 



(5.8) 



where wt{p) is the transverse mode dispersion relation and is given in Refs.[14, 21], together 
with pI'^^. Taking Kq > m + mg, without loss of generality, we obtain 



/P 



ole 



8nnB{Ko 
■nBiwtipo)) 



^pr{wt{po),Po)nB{Eo) 

-C/O 



1 



dE/dp"^ + dwi/dp"^ 



(5.9) 



P=PO 



with po given by the 5-constraint: 



Wt{po) = Kq - Eq = Kq - \pI + m 



1^ _L 'rr)2 



(5.10) 



In order to perform the interesting limit m — > and Kq — > m^, we notice that the solution 
of (5.10) vanishes for to = as 

Pq = Kq -mg + 0{pl), 



usmg 



wt{po] 



Po 



^ 5 nig 



(5.11) 
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Working out the derivative with respect to m^ and taking m ^ 0, we find that the most 
singular term in the hmit po ^ comes from the derivative 



This allows to derive in the appropriate limit: 



ole - — 1 ^ 1 

A"°'K,0)^^^^^-^^^^— — ^. (5.12) 

In order to conclude about /i, we should now take care about the contribution coming 
from the branch cut of the effective gluon propagator associated to Landau damping. In 
this case, from first inspection, it seems safe to work directly on-shell and to calculate 

/f-K,0)=Um(-^)|-=. (6.13) 

Going back to Eq.(5.7), we end up by writing if^^"^ as the sum of two terms 1°^ and /^, 



-d^J j,±:^ P'dp^nsim, - E)pf-{m, - E^p\ 
/^oc hm f-/^) r p2^pM^(l + nB(i?-m,))pf-(i?-m„p), 



(5.14) 



neglecting irrelevant numerical factors and with p^^*^*^, the transverse spectral density as- 
sociated to Landau damping, given in [21]. 

When taking the deriative -^^ each term yields two contributions: 

- the first one is given by the derivative of the integration bounds. In this case, the only 
source of potential singular behaviour is associated, in both /" and /^, to the value of 
the integrand at p = E = nig. Since lim pf^^'^{w,p) ■ nsiw) is finite, no singularity is 
encountered. 

- the second one is obtained when taking the derivative of the integrand. Writing 
dm^ ~ ^'§E ^'^'^ integrating by parts, we end up as above with a finite contribution. 

As a consequence, the contribution of Ii{k) to Im A*nt(/c) is, in the appropriate limit - 
neglecting numerical factors, 

4 



Im A*n/'^) oc ig-'N— .^^ ^^ .„ . (5.15) 

rUg [Ko - rUgY 
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With 



kT = ml- Kq^ 2mg{mg - Kq) 



this leads to 

Im A*Ut^^'^ oc Cg^NTrUg . (5.16) 

(ii) Calculation of the longitudinal integral 

In order to complete the argument, we need to calculate l2- This in fact, using 
Po = p'^ ~ P^j amounts to calculating two integrals: 

/(IT) T) 
-^^^^2 , (5-17) 

soft 

and 

hB = Im 1^4^. (5.18) 

J Q P Pl 

soft 

These two integrals may be calculated, following the same line as above, from the basic 
integral 

I(k) = Im / ^5 ^^5 7^^-o J (5.19) 

J q + TO p + m'^ pj^ 

for k = 0. This integral, analogous to a vertex integral, is evaluated with the technique of 
Ref. [8], already used here which starts by introducing "mixed" propagators for which an 
integral representation in terms of spectral densities is given. To be specific we introduce 
the mixed propagator A(r, q) as: 

A{qo, q)= f rfre^^°^A(r, q) = ^-^ , (5.20) 

Jo Q^ + m^ 

and similarly for A{po,'p) and A^^p) = 4r. 

Pl 

This allows to write I{k) as 

I{k)=Im J-^J dnj rfr2e^^i'="A(ri-r2,i?p)A(ri,£;,)Ai(r2,p) , (5.21) 



with Ep = a/p^ + m^. Eg = \/p^ + m'^, and in the following p^ = p 
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The spectral representation which is then used in order to calculate the imaginary 
part, 



+ CO 



A{t,E)= [ dw e-'"^p{w,E){l + nB{w)) (5.22) 

— oo 

is only valid when < r < /3 and therefore A(ri — r2, Ep) should be defined to be periodic 
with period 13 [8]. This leads to (after analytic continuation: z/cq — ^ Kq + ie) 

I{Ko, k = 0) = 

+ 00 +CXD +00 



{2'kY J J J W — W2 

— oo — oo — oo 

■ nB{wi){nB{w2)d{Ko - wi ~ W2) - nB{w)d{Ko - wi - w)}. 

The calculation then follows closely the steps described for /i. Let us first consider /2a and 
focus on the pole contribution in the spectral function for the longitudinal mode pl{w,p): 

pI^'%w,p)=pT{w,p)5{w-wl{p)) = -^ ji'^'~/'\ siw-WL{p)). (5.24) 

p^ Snig — w^ + p"^ 

with wl{p) the dispersion relation given in Refs. [14, 21]. 

The only interesting configuration in the limit Kq -^ nig is due to the term propor- 
tional to d{Ko — Eg — wl{p)) in complete analogy with the case of the transverse integral. 

We find 

00 00 

^^(^.o) - .rr^ (-^] f'-^P' flu. 



8nnB{Ko) \ dm? J J Eq j 

^0 (5.25) 

■ ^^Q^^nB{Eq)nB{w2)5{K^ - E^ - w^) . 

It turns out that there is a complete correspondence with the calculation of /i, i.e. replac- 
ing: 



Pt {wt[Po),Po)^ 2 ^TT^ ' (5-26) 



Po PTi^LiPo),P o) 

pI-wKpo) 

which both equal 2^ when pq -^ 0. We obtain therefore the same singular value for /2a 
and /i. For /2b, the additional differentiation with respect to m''^ amounts to replacing 
the above correspondence by: 

res/ / X X PtPTi'WL{Po),Po) /^^„^ 

Pt iMPoJ^Po) -^ -Y^ 2/ ^^2 ' ^-27 

{po-wi{po)y 
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which in the hmit po ^ yields l2b — —^Ii- The integral l2b is therefore negligible in 
the relevant on-shell limit. 

Going back to (5.2), we find that the gauge dependence of the transverse structure 
function is indeed non vanishing on the mass-shell 

Im A*Ut{mg,6) = -ig^Nk^ {h{mg,Q) + Ia{mg,Q)) oc g'^NTrUg , (5.28) 

which implies a gauge dependence for the damping rate (cf. Eq.(3.21)): 

S^^ oc ^Jll^Jk oc ^g^NT , (5.29) 

rUg 

for the transverse gluon excitation at rest. 
Coulomb gauges 

In this case one may vary the effective propagator (Eq.(4.4)) as 

*D^'^ ^ ^Df"" +Cc^^ (5.30) 



r2 ,2 
can be written as in (4.6) with A/3 = and Act = ^^ \ so that the gauge variation 



around the strict Coulomb gauge (^c = 0) . The gauge dependent part of the propagator 
can be written as in (4.6) with A/3 = and Act = 
Im A*Ilt is obtained from (5.2) with trivial changes: 

Im A*Ut - -Ccg^N Im 



dp I p"^ {k ■ p)'^ 

soft ^ V ~1^^ Pp^ 



-^TJ „S|l-^ + 






4 ^JsoitQ^P^ \ Pp^ 



/j.2p*2 

The only potentially dangerous contribution is proportional to the integral 

dp 1 



Q^PlP^ • 

soft 

It does not depend, however, on ko and therefore this yields - after continuation - a van- 
ishing imaginary part. 

We thus conclude that the problem of gauge dependence arises in the general class of 
covariant gauges but not in Coulomb gauges. The implications of our results are discussed 
in the next Section. 
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6. DISCUSSION 

In the previous Section we have shown by exphcit calculations that the gluonic damp- 
ing rate for the excitation at rest remains gauge parameter dependent in the class of 
covariant gauges. This extends the result previously derived for the fermionic rate [11] in 
the same framework of the resummed effective perturbative expansion [8] . 

In the course of the calculation in Section 5 it is stated that the integrals in Eq.(5.2) 
have no imaginary part on the mass-shell Kq = nig for all values of the mass parameter 
m > 0. This is most easily seen from the constraint (5.10) which cannot be fulfilled when 
Wf = nig] a solution is only possible on the mass shell for ni = 0. One might therefore 
argue that continuity in m requires the expression of Eq.(5.2) to vanish on shell for m = 
as well; indeed this would be the case if the order of limits (cf. Eq.(5.5)) were interchanged 
to lim^^o liniKo^m • However, the parameter ni was introduced as a technical device 
to treat the double pole 1/(13'^)^ in the gauge dependent term in Eq.(5.4), so the order of 
limits is fixed as in Section 5. Nevertheless this observation suggests that the observed 
gauge dependence may be intimately related to infrared divergences, despite the fact that 
the damping rate itself is definitely infrared finite. 

The discussion may be clarified by introducing an infrared regulator, as advocated by 
Rebhan [22] (see also refs.[23] and [24]). In order to demonstrate the problem, we consider 
in detail the integral (cf. Eq.(5.6)): 

Hko.k = 0)= /-^, (6.1) 

J 1 FT 

which is an intermediate step in deriving a contribution to the self energy of the same form 
as the first term in the second line of (5.2), namely 



/(fco,fc = o)~4 / rrf^- (6-2) 



4 / ^ 

(q4) pI 

The discussion for these integrals may be generalized to the more complicated cases cor- 
responding to the integrals in Eqs.(5.3), (5.17) and (5.18), and even to A*Ilt (Eq.(5.2)) 
itself, before taking the imaginary parts. However, in order to simplify the analysis, the 
simpler dispersions [9] p'^ — Pq+P ^ + "^3, and p | = n'^p'^ are used, instead of the physical 
dispersion relations (Eqs.(3.9, 3.10)). 

We start in the imaginary time- framework, /cq = 2nniT, and later the continuation 
iko -^ Ko+ie is performed. In order to regularize possible infrared singularities dimensional 
regularization is applied with analytically continued spatial dimensions n = 3 + 2e [23] : 

"p-T j: j^^. (6.3) 

P()=27rmT 
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In the high temperature hmit the leading term of (6.^ 



is: 



T /(A;2 + m2)2 r^ p2. 

/(fco,fc = 0)oc^= X \'' / dp ^ 



A^V ^l Jo V + i^iS^ 



T 1 1 

oc^ r(-+e)r(- 






4fc„ 



(6.4) 



Irrelevant numerical factors are neglected. As long as the energy /cq is discrete and euclidean 
there are no singularities present and we are allowed to take e = 0. In the next step the 
analytic continuation to the Minkowski energy Kq is performed, and the integral becomes 
imaginary: 

Imi{Ko,k = Qi)(x-— , (6.5) 

and independent oinig. One may note that at non-zero temperature Im I is also defined 
for Kq < nig, in contrast to the T = contribution, which has only support for Kq > nig. 

The integral I{ko,k = 0) of (6.2) contains a term ~ k'j^ T L dp p^'^~^, which is 
infrared singular in three spatial dimensions. Therefore a non- vanishing regulator e seems 
appropriate, at least for the real part of the integral (after analytic continuation) . 

However, it is remarkable that indeed these infrared singularities cancel in A*nt - at 
least in the described approximation: A*nt has precisely the same structure as exemplified 
by the integral in (6.4) ! 

The imaginary parts of the integrals (6.1) and (6.2) are smooth functions in terms of 
the continued external energy; for e = and at the value of Kq = nig we find e.g. for / 

from (6.4), 

T 

Ini I(nin, k = 0) (X — , (6.6) 

nig 

in agreement with the result for /p°^'= in Eq.(5.9) for Wt = trig and m = 0. I.e. a non- 
vanishing answer is obtained with the following prescription: first the necessary integra- 
tions are performed for discrete external euclidean energies, where an infrared regulator 
may be introduced, which is then removed before the analytic continuation to external 
Minkowski energies is done. The discontinuity of the integrals has to be evaluated when 
staying off-shell, and finally the on-shell limit is taken. This prescription is equivalent to 
the calculation of the imaginary part of the integrals as described in Section 5; it appears 
to be the appropriate procedure when dealing with the effective euclidean theory at finite 
temperature as it is implemented by Braaten and Pisarski [8]. However, it leads to a 
dependence on the choice of the covariant gauge parameter. 
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The resolution of the gauge dependence problem advocated in Refs. [22-24] follows the 
steps of taking the on mass-shell limit just after the continuation, but keeping the infrared 
regulator e different from zero. As emphasized by Rebhan[22], the correct treatment of 
mass-shell singularities at zero temperature requires an infrared regulator in the framework 
of perturbation theory [18]. For the example above it is obvious (cf. Eq.(6.4)) that 

I{Ko=mg,k = 0) = 0, (6.7) 

for e > 0, so there is no imaginary part at Kq = nig, and gauge independence is restored! 
Although this prescription does appear to fix the problem of gauge dependence at this 
order, it is perhaps appropriate to make the following comments: 

1. At zero temperature, the infrared regulator is required to make the wave function 
renormalization well defined. In order to evaluate the pole position, the regulator is 
not required. In the present case, we are only interested in the pole position and in the 
effective (truncated) approximation under consideration for calculating damping rates 
the necessity of the infrared regulator is not at all obvious: only the imaginary parts of 
the integrals in (5.2) are physically relevant (since higher loops can contribute to order 
g^T to the real part of the self energy) and these imaginary parts are not infrared 
divergent. (We repeat that even the real part of A*!!^ is free of infrared singularities 
in the approximation investigated above.) Therefore the limit e ^ can be taken in 
(6.4) without encountering a singularity. Using the prescription in [22,23], however, a 
discontinuous behaviour of the function Im I{Kq, k = 0) (and therefore of A*!!^) is 
enforced: it vanishes at Kq = nig, whereas for Kq ^ nig it behaves as T/Kq (in the 
limit that e = 0). In this way, the value of the damping rate, a supposedly physical 
quantity, is affected by the presence, or absence of such a regulator. Clearly, the 
question as to whether or not the regulator is necessary in the expressions considered 
is considerably more subtle than at zero temperature. 

2. As long as AX in the gauge dependence identity (1.1) is well-behaved on mass-shell, 
gauge dependence is guaranteed algebraically order by order in any self-consistent 
perturbative expansion. In the presence of the mass shell singularities discussed above, 
the gauge dependent terms in the lowest order damping rate need to be "regulated 
away". They do not vanish algebraically, but instead rely on the presence of the 
regulator to control the value of the integrals in question. This mechanism for ensuring 
gauge independence is very different in this regard from the algebraic proof, and it is 
not clear whether it will work at higher orders. (This remains an important question. 
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in our opinion, even if such higher order contributions are difficult, if not impossible, 
to calculate in practice.) 
3. For gauge independence to hold, the infrared regulator must preserve the Ward iden- 
tities of the theory. Dimensional regularization[23] does so trivially, but results in a 
somewhat strange analytic behaviour for the off-shell propagator due to terms of the 
form {Kq — m^gY. While a cut-off avoids this particular problem[22], verifying the 
Ward identities is somewhat more problematic. 

To summarize, one appears to be left with the following choice: one can retain an 
infrared regulator throughout the calculation, with potential consequences as outline above, 
or one may use the argument of continuity (this time in Kq) to follow the prescription 
advocated in Section 5 for calculating the damping rates in the effective theory of Ref.[8]. 
However, if the latter choice is made, it remains an open problem to determine whether 
or not there exist missing terms [25] in the resummation which restore (covariant) gauge 
independence of the damping rate. 
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